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ABSTRACT 

If accretion disc emission results from turbulent dissipation, then axisymmetric accretion the- 
ory must be used as a mean field theory: turbulent flows are at most axisymmetric only when 
suitably averaged. Spectral predictions therefore have an intrinsic imprecision that must be 
quantified to interpret the variability exhibited by a source observed at different epochs. We 
quantify contributions to the stochastic imprecision that come from azimuthal and radial aver- 
aging and show that the imprecision is minimized for a particular choice of radial averaging, 
which in turn, corresponds to an optimal spectral resolution of a telescope for a spatially un- 
resolved source. If the optimal spectral resolution is less than that of the telescope then the 
data can be binned to compare to the theoretical prediction of minimum imprecision. Little 
stochastic variability is predicted at radii much larger than that at which the dominant eddy 
turnover time (~ orbit time) exceeds the time interval between observations; the epochs would 
then be sampling the same member of the stochastic ensemble. We discuss the application of 
these principles to protoplanetary discs for which there is presently a paucity of multi-epoch 
data but for which such data acquisition projects are underway. 

Key words: accretion, accretion discs; (stars:) planetary systems: protoplanetary discs; tur- 
bulence 



1 INTRODUCTION 

Accretion discs are ubiquitous in astrophysics, often forming 
around stars and compact objects where angular momentum would 
inhibit direct gravitational infall of the plasma onto the central ob- 
ject. A long standing theme of research in accretion theory has been 
to understand how how the discs transport angular momentum and 
remain quasi-steady accretors. This transport likely involves differ- 
ent combination of local and global mechanisms depending on the 
circumstances. However, there is presently an intellectual gap be- 
tween the study of the fundamental theory of angular momentum 
transport in accretion discs via numerical simulations and practical 
accretion disc models that can be used to compare with observed 



spectra. For reviews that sample a range of perspectives see Frank 
[etaLl p002]> ; |B albus and Hawley | < f2003l > ; |Hartmanii| ( [2009) ; Spruit 
( |201ot ; |Blackman| ( |2010j i~ 

The most commonly used practical accretion disc model is 
based on local viscous transport of angular momentum in axisym- 
metric discs and hides the unknown physics of transport into an 
effective local "turbulent viscosity" ( |Shakura and Sunyaev||1973| 
Pringle 1981 ; Frank et al. 20021. Because microphysical molecular 
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viscosity is typically much too small to drive the observed accre- 
tion rates given constraints on surface densities, some kind of en- 
hanced transport is needed. While these mechanisms routinely in- 
volve some form of turbulent motion (often with a magnetic origin), 
the standard practical model equations describing accretion discs 
do not specify or accommodate the nuanced physics of transport, 
lumping the processes into a single dimensionless transport coeffi- 
cient. The assumption that the discs are axisymmetric can then at 
best apply in a spatially or temporally averaged mean-field sense 
since turbulence necessarily violates axisymmetry locally and tur- 
bulent transport is a stochastic process. 

In order to properly compare standard axisymmetric accre- 
tion theory with observations, the limitations of the theory must be 
quantified. The theory has a limited precision which implies that 
agreement between models and observation can only be approxi- 
mate. Disagreement within the precision error may be consistent 
with a stochastic nature of the discs. A key quantity, and the focus 
of this paper, is the relative precision error (RPEQin the lumi- 
nosity at a given frequency. Blackman] {1998| l discussed the RPE 



1 Compare the statements 'babies sleep between one and twenty three 
hours per day' and 'babies sleep 8± 1 h per day.' The first statement has 
a large RPE, but is almost certain to be accurate. The second statement has 
a low RPE, but may well be incorrect. However, its low RPE makes the 
second statement testable 
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for accretion discs in the context of Advection Dominated Accre- 
tion Flows (ADAFs) and thin discs in active galactic nuclei (AGN) 
and X-ray binarys (XRBs) in which the observation times are often 
longer than the orbit times at a given radius. In contrast, protoplan- 
etary discs often involve small observational exposure times com- 
pared to disc orbit times. In this context, we revisit quantifying the 
precision of standard accretion disc theory and make predictions 
for multi-epoch observations. The predictions can help to diagnose 
whether a systematic variability may also be present. 

In section 2, we compute the azimuthal and radial contribu- 
tions to the luminosity RPE and combine them for a spatially un- 
resolved source. In section 3 we discuss how to minimize the RPE 
and the role of telescope spectral resolving power. In section 4 we 
consider how to apply the results to the specific case of LLRL 31, 
one of the few protoplanetary discs in which there are multi-epoch 
observations. Although we find that the observed variability in that 
source is systematic and not stochastic, the application exemplifies 
how to use the results herein. We also discuss further implications 
for future multi-epoch observations, of protoplanetary discs, and 
conclude in section 5. 



CALCULATION OF RPE IN LUMINOSITY FROM 
TURBULENT DISSIPATION 



We assume a standard a-disc prescription (Sha kura and Sunyaev| 
|1973[ > for the turbulent viscosity.The characteristic turbulent eddy 
scales have characteristic velocities v csl and sizes l csl so that the ef- 
fective viscosity satisfies 



v ~ ac s h ~ w c d^e, 



(1) 



where c s is the disc sound speed, h is the scale height, and is 
the orbital speed. We take < a < 1 to be a constant over the 
whole disc. The last similarity in {T} follows from the assumption 
that the eddy turnover time r e£ j ~ i consistent with simulations 
of the magnetorotational instability (MRI) (e.g. Balb us and Haw-| 
|ley| ( |2003| ). For a thin disc in hydrostatic eq uilibrium, Eq. (ffl> also 
implies that « ec j ~ a 



and l e d ~ a^^h (Blackman 



q. fj); 

1998) . 



In what follows, we focus on cases for which the required ob- 
servational exposure time is much less than the eddy turnover time 
at all radii. We can then regard a measured spectrum as an instanta- 
neous snapshot. This applies to Spitzer observations of protoplan- 
etary discs. The discs range in size from a few tenths of an AU to 
several tens of AU so that t ed ranges from days (in the inner disc) 
to years (in the outer disc). In contrast, a typical Spitzer exposure 
takes only a few minutes. 

Since the sources of interest are spatially unresolved, we cal- 
culate the contributions to the RPE by averaging in the azimuthal 
and in the radial directions for the axisymmetric theory. If the disc 
is sufficiently optically thick then the photosphere occurs within 
the first layer of energy containing eddies. An optically thin disc 
requires additionally averaging over the third dimension of eddies. 



2.1 Azimuthal Contribution 

Consider an annulus of width equal to one eddy scale centered at a 
radius r within the disc. The number of eddies producing observ- 
able emission in this annulus is given by 



2nr 



1 + 



h 

led 



(2) 



where t v is the optical depth. In the limit that t v — > 0, eddies 
over the entire thickness of the disc contribute to the emission. For 
r„ — > oo only the top layer of eddies contribute to the emission. 

Over its lifetime, the brightness of each dominant eddy will 
rise to some maximum and then fall as it cascades to small scales. 
Let us assume that the rate of brightening and dimming are constant 
so that the eddy brightens to its peak luminosity Lo at t e d/2 and 
then falls to zero at t = to- The probably distribution function for 
the luminosity during the rise and fall is then 



p(L e )dL e 



1 



Then the mean luminosity for each eddy is (L e ) 
mean squared luminosity is 



(Li) 



+ 



The variance is then 



\L e ) = (Ll) - <L e > 2 = 



Mi 
12 



(3) 

Lq /2 and the 
■■ io/3. (4) 

(5) 



satisfies 



y(L) 



If the value of the mean luminosity of each eddy is normally 
distributed about Lo/2 then the standard deviation of the luminos- 
ity per eddy when averaged over all N e d eddies in a given annulus 
Using L v to indicate the luminosity at frequency 



v measured from the annulus whose peak temperature corresponds 
to that frequency, the fractional luminosity variation of the entire 
annulus of N e d eddies is then 



N ed AL e 
N ed (L e ) 



a(L e )/VN 
L /2 



(6) 



where we have used Eq. {2|, the relations below Eq. and the 
subscript az to indicate the azimuthal contribution to the variation. 

r -i 1 1 

In the optically thick limit Eq. 1 6 1 gives y ^ J = ^= ( £ ) 2 . 



2.2 Radial Contribution 

In determining the radial contribution to the luminosity RPE we 
consider two independent contributions. The first comes from a ra- 
dial smoothing implicit to the mean field theory. The second comes 
from turbulent velocity fluctuations that add macroscopic fluctua- 
tions about the mean orbital velocity at each radius. We address 
each in turn. 

In formulating a mean field theory for accretion discs, one can 
in principle simply take vertical and azimuthal averages to obtain 
a theory that depends only on radius without any radial smoothing. 
However in the commonly employed framework which invokes a 
turbulent viscosity (Shakura & Sunyaev 1973), the scale of the tur- 
bulence is assumed to be smaller than the scale of the mean quan- 
tities. The theory is not really resolved on spatial scales below an 
eddy scale so we therefore include a radial smoothing scale, £, cho- 
sen such that Z ed < £ < raise . In addition, telescopes have a finite 
spectral resolution and a given frequency in practice corresponds 
to a radius range in the disc, not a precise radius for a spatially un- 
resolved source. We therefore invoke a radial smoothing of mean 
quantities of the form 

X(r) « (X(r, t)) t = J ^ X(r + A, t)d\, (7) 

where X is an arbitrary quantity to be averaged and X is its mean. 



Imprecision of Accretion Theory and Protoplanetary Discs 3 



The similarity follows from the assumption that eddy statistics are 
steady in time. Blackman (1998) discussed an additional temporal 
average but as mentioned above, our present interest is for cases in 
which the required snapshot observation time is much less than an 
eddy turnover time at all radii. 

The radial contribution to the luminosity variation can be writ- 
ten as 



Disc Model 



" AL„" 
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rad 


L v dr 
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r 



(8) 



We must now calculate the contributions to the right hand side. 

We write the total imprecision in radius as a quadratic sum of 
the smoothing and velocity fluctuation contributions. That is, 

(Ar) r 2 s 



(Ar) 



+ 



(Ml 



(9) 



where first term comes from the blurring induced by the radial 
smoothing scale £ and the second term comes from the macroscopic 
turbulent velocity fluctuations. 

For the radial smoothing term of Eq. {9}, we use the fact that 

r ^ <r <r+ -^j- (10) 

2NS 2iV rs 2 

are indistinguishable once £ is chosen, so that 

(Ar) rs _ £ 



2rNi 



where 



7V rs = 1 + 
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(12) 



and measures the effective number of averaging scales per f . The 
role of t„ the same as that for Eq. ([2j. In the limit t a ba « ted and 
r„ >> 1, Nfi = 1. For the velocity fluctuation contribution to Eq. 
l|9j, we note that the dominant contribution to the local velocity at 
a given radius is oc r _1//2 from Keplerian rotation. Thus 
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where we have again used the standard result for error on the mean 
as applied to Nfi eddies each of whose mean velocity is normally 
distributed with fluctuations of order v e d- The number of eddies 
over which the radial averaging is performed is given by 



N„ = 



- 1 



+ 



^obs 



(14) 



For r„ >> 1 and t bs/ted ~ 0, this gives Nfi = £/l e d- In this 
limit, adding {TT} and |T3j in quadrature, we obtain 



(Ar) 



,2. a 2 

v 2 £ 4r 2 ' 



(15) 



If each ring of material radiates around a specific blackbody 
temperature T, the luminosity of an annulus of width l e d is 



where 



L v = 2jrrfl?* (T(r)) l e 



B p J dk (T) = AT 3 



(16) 



(17) 



is the peak value of the Planck function B v and A is a constant that 
comes from evaluating 



Bu{T) 



2k A 



h 2 c 2 e x 



(18) 



Viscously Heated, Optically Thick 9/10 
Viscously Heated, Optically Thin 3/5 
Irradiated by central star 1/3 



Table 1. Values for p in Eq. j20| for various disc models I Ed gar et al.|2007) 



at x = hv/kT ~ 2.82, corresponding to the Wien displacement 
law. Then 



dL v 
dr 

If we assume that 
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then combining with the scalings below Eq. (1) gives l e d = 
a 1/2 h oc c s /Q oc r 5 (3_p) . In combination with Eqs |s}, |l6 
and|l9), we then obtain 
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|Edgar et al.|p007| > calculated p for a variety of disc models, and 
their values are given in Table [T] Eq. \2\ \ is strictly applicable for 
the Tv » 1 regime so here p — 9/10. 



3 TOTAL RPE AND ROLE OF RESOLVING POWER 

Eqs. ([6| and \2l\ give the respective azimuthal and radial contribu- 
tions to the RPE of the predicted luminosity in terms of disc model 
properties. They must be added in quadrature to get the total lumi- 
nosity RPE. The choice of £ can be taken to minimize the RPE but 
whether this minimum can be compared with data depends on the 
resolving power of the instruments. We explain these points further 
below. 



3.1 Minimizing the Intrinsic RPE 

First we determine the minimum of Eq. (15) , which applies in the 
optically thick r v » 1 and t t s = limit. Differentiating and 
setting equal to zero gives 



Sopt 
r 



(22) 



The optically thin case gives the same result because the optically 
thin corrections to N rs and Nfi just produce an overall multiplica- 
tive factor of l e d/h to both terms in Eq. 1 15 l and thus do not change 
its minimum. 

Choosing £ = £ opt will minimize the RPE of the theory. In- 
creasing £ would average over more turbulent fluctuations, but at 
the price of poorer spatial precision. Reducing £ would increase 
the random error due to a smaller number of eddies. If £ opl < i ed , 
we must use £ = l ei instead. 



Substituting Eq. (22) into Eq. {T5}, gives 



Ar 



(23) 



where we have used v^l^ — ac s h. Note that £ opt of Eq. \22\ can 
only be used to compare with observations if the observations are 
resolved on scales £ ^ £ op t since only then can the data be binned 
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accordingly. We will come back to this point. First we calculate the 
RPE for the luminosity using £ = £ pt- 

Substituting l e d = a 1//2 fa and v e d 
and in turn substituting the result into Eq. l |2 1 [ i gives 

AL V 



a 1/2 c s into Eq. 1 23 i 
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(24) 



Combining results from Eqs. |6]l and l |24) by adding the radial and 
azimuthal contributions in quadrature,we have for the net minimum 
RPE in the r v » 1 limit 

1/2 

+ 



AL V 



(25) 



where in the optically thick limit (Edgar et al. 2007 p = 9/10 and 
h/r is given by 



- = 1.8 x 10~ 2 
r 



M, 

lO-^Qj/r" 1 



M, 
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(26) 



0.04 



0.03 



0.02 - 



0.01 - 



o.oo L 




500 1000 1500 

Resolving Power, R 



2000 



Figure 1. RPE in luminosity from turbulent dissipation for 
optically thin case at r = 10AU for the parameters of 
LLRL 31 as a function of resolving power R e f for a = 
0.1(Thick lined Blue), 0.01(thin lined Red), O.OOl(Dashed). The higher 
a has the most prominent minimum, which does not correspond to the 
maximum resolving power but an intermediate value. 



3.2 Expressing Total RPE as function of Resolving Power 

For a telescope, we can only choose £ to be as small as the spatial 
or spectral resolution allows. If the smallest £ that the telescope 
can resolve is smaller than £ op t, the data can be binned to create an 
effective £ = £ opt that minimizes the RPE. 

We can estimate the radial averaging scale from the spectral 
resolving power of the telescope. The resolving power is given by 



R : 



v 

Av' 



(27) 



From the Wien displacement law, we know that a blackbody spec- 
trum has i^peak oc Tpeak, so the effective resolving power is then 
R b} « T/AT. From Eq. J201, we then obtain 



Substituting into Eq. (15) , we find 
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(29) 



Using this in |2TJ we obtain 

AL V 
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5 _ 7p 
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(30) 

The azimuthal contribution Eq. l|6j, does not depend on the fre- 
quency explicitly so combining it in quadrature with Eq. < |30| > gives 
the total effective RPE 



data are not more coarsely binned by hand, the data could represent 
a radial resolution less than the eddy scale of the theory. The a-disc 
model based on turbulence becomes ill-defined below such scales. 
Coarse graining of the data has to be considered carefully before 
sensibly comparing with a mean field theory. 



3.3 Optically thin case 

Finally, note that for t v < < 1 Eqs. j2]l ( |12) and l |14) take a different 
values compared to the t„ >> 1 case. The overall calculation of 
the RPE follows that which led to |3TJ with the analogous result 
being 



7p|2 
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(32) 



where in this r v « 1 limit (Edgar et al. 2007 1, p = 3/5 and 
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(33) 



To obtain minimum (optimal) value of (32) (the analogue of (|25)), 
we substiute in l |32) the value R e f = R opt = r/p£ opt and use 



Copt jr from Eq. 1 22 1 which applies for both the optically thin and 
thick cases, 



AL„ 



1£ 

2 



Aa^ 2 ^ypR Bf + 



4p 2 R 2 



+ 



1/2 
67T 
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(31) 

It might seem surprising that increasing R in Eq. ( |30) can give less 
precise results but recall that for a fixed p and r there is a one-to-one 
correspondence between R e f and £. But as discussed above, there 
is a value £ = ^ opt that minimizes the RPE. There is therefore a 
corresponding value of R e f that minimizes the RPE and the two 
turn out to be related by £ op t = pfl r - . Figure 1 illustrates these 
basic principles in plots of Eq. ( |31) for different values of a. Note 
that the minima are less pronounced as a is decreased. 

Note that if R of the telescope is sufficiently large and the 



.4 IMPLICATIONS FOR MULTI-EPOCH OBSERVATIONS 

A basic application of the RPE described above arises for multi- 
epoch observations of a single object when the time scale for a 
snapshot observation is shorter than the eddy time at any radius. 
Emission from radii at which the time scale between epochs of ob- 
servation exceeds the eddy turnover time would then be expected 
to vary stochastically between epochs. In contrast, at radii larger 
than the that at which the eddy turnover time matches the epoch 
interval time, successive observations sample the same member of 
the statistical ensemble and no stochastic variability corresponding 
to the the local energy-dominating eddies would be expected. 



Imprecision of Accretion Theory and Protoplanetary Discs 5 



0.015 - 



0.010 



— J 
< 



0.005 - 



0.000 



4 6 
Radius, r (in au) 



10 



Figure 2. RPE in luminosity from turbulent dissipa- 
tion vs. radius r for the optically thin case, using a = 
0.1(thick lined Blue), 0.01(thin lined Red), O.OOl(Dashed). From left to 
right the three vertical dashed lines correspond to three specific radii whose 
orbit period equals different inter-epoch observation periods for LRLL 3 1 
(Muzerolle et al.|2009) : 1 week (r = 0.087 AU) and 4 months(r = 0.58 
AU), and ~ 25 months (r = 2 AU). For the respective inter-epoch periods, 
these vertical lines mark the upper radial limit of observational validity 
of the curves, and above which the expected variability will be sharply 
reduced. For example, if the inter-epoch interval is 25 months, then the 
RPE would predict the regime of validity for the curves to be r < 2au. 
Above r > 2 AU the prediction would be that AL„/L„ falls near zero 
because the epochs would sample data from the same statistical ensemble. 



4.1 Lessons from the specific case of LRLL 31 

Presently, multi-epoch data are sparse. Programs such as YSOVAR 
(Stauffer et al. 2010) will eventually produce more data. Neverthe- 
less, to focus the practical application on a on a specific object, we 
consider IR Spitzer multi epoch observations of the source LRLL 
3 1 in the star-forming region IC 348 l |Muzerolle et al.|2009) . For 
reasons discussed below, this object is most likely not showing its 
primary variability from stochastic effects we discuss but our anal- 
ysis to reach this conclusion below illustrates how to use the con- 
cepts of the present paper to arrive at that conclusion. 

The stellar mass of LRLL 31 is M = 1.8Mq and the es- 
timated accretion rate is M = 1.5 x lO~ 8 Af0yr _1 (from Pa/3). 
For the low spectral resolution mode used to observe this source in 
the 5-38 micron range (Muzerolle et al. 2009), the Spitzer spectral 
resolving power is R ~ 200. The observation times are minutes, 
well below the smallest relevant orbit time and consistent with set- 
ting t bs ~ in \12\ and l |14) and out assumption throughout the 
previous sections. 

Figure 2 shows Eq. < |30| > plotted as a function of radius for 
three different values of a and a value of R c f = 200 for the opti- 
cally thin case. The vertical lines show, for three separate choices 
of the inter-epoch time, the radii above which there would be no 
variability expected because the inter-epoch time is less than the 
eddy turnover time above these radii. 

Tables 2 and 3 exemplify how use of Ref = R, the resolving 
power for the telescope, compares to the use of R e f = R op t, the 
value that minimizes the luminosity RPE, for different choices of r 
and a. For these tables, the epoch time is taken to exceed the eddy 
turnover time at all radii. In reality, as emphasized by the vertical 
lines in Figure 2, there would be negligible expected variability at 
radii above which the orbit time (~ eddy turnover time) exceeds 
a chosen inter-epoch observation interval. For cases in which the 
telescope resolution (column 6) exceeds the optimal value (column 



4) the data can in principle be binned such that the binned data can 
be compared to the optimal a-disc model prediction. The tables 
highlight two trends in the Tables: (i) the larger the value of a the 
larger the predicted variability between epochs at a given radius 
and (ii) the larger the radius for a fixed a the larger the predicted 
variability. 

Muzerolle (2009) argue that the disc of LLRL 31, at r < 10 
AU, is optically thin and substantially cleared out (perhaps via the 
action of planet) below 10AU but optically thick for a least some 
range of larger radii. Thus for this object, the values in the Table 
2 would apply for r < 10AU and the values of table 3 may apply 
for some range of larger radii above which the disc again becomes 
optically thin. The value at which the latter occurs requires more 
detailed modeling. 

There is also a transition radius above which viscous dissi- 
pation would be dominated by irradiation from the central star 
( |Sicilia-Aguilar et al.|20"06l|Edgar et al.|2007fr , namely 



0.3AU 



I 



(34) 



10- s M Q yr- 1 / 

For r > r c the photospheric temperature induced by viscous dissi- 
pation falls below that associated from surface dust illumination by 
the central star. The values of p would then flatten to p — 1/3 
(see Table 1) and irradiation would dominate at all larger radii. 
Disc emission resulting from stellar irradiation is less influenced 
by the stochasticity of turbulence than emission via turbulent dis- 
sipation. Since the local blackbody emission varies with T 3 , we 
would expect a much reduced stochastic variability in regions of 
the disc dominated by emission from reprocessed starlight. This 
implies that even for epoch intervals longer than the orbit periods 
at r > r c that the stochastic variability discussed herein would be 
reduced compared to regions of r < r c It is important to emphasize 
that the RPEs calculated in Tables 2 and 3 is a stochastic variabil- 
ity corresponding to that of mean field a-disc resulting from the 
disc luminosity contribution by stochastic process such as turbu- 
lent viscous dissipation. As applied to the specific object LLRL 3 1 
therefore, these would only amount to a variability in a subdomi- 
nant contribution to the total emission above r c . 

An absence of observed stochastic fluctuations for r > r c 
would support the expectation that the primary emission in those 
regions is not the result of turbulent dissipation. For r < r c the ab- 
sence of stochastic variability would be expected only if the inter- 
epoch time scale is shorter than the characteristic eddy turnover 
time at the radius dominating the radiation at the particular fre- 
quency measured. Note also that the amplitude of the predicted 
stochastic variability would not depend on the epoch interval as 
long as the epoch interval is larger than the orbit time at the ra- 
dius producing the corresponding emission. Other systematic non- 
axisymmetry or variable accretion rates would add to the variabil- 
ity and could be distinguished form stochastic effects given enough 
multi-epoch observations 

For the specific case of LRLL 31, the largest time scale be- 
tween epochs is of order 4 months which means that according 
to Fig. 2, the stochastic variability would be expected only for 
r < 0.58 AU. This is well within the dust depleted inner region 
of this transitional object where the disc is optically thin and where 
the emission would come from dissipation internal to the disc rather 
than reprocessed starlight. Thus there would be some stochastic 
contribution to the measured variability between epochs predicted 
by the RPE of the mean field models. However, the double-digit 
percentages of the observed variability reported by Muzerolle et al. 
( 2009 1 in LRLL 3 1 are larger than the values predicted in our Ta- 
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bles that arise from stochastic variability alone In addition, the 
presence of a single pivot wavelength of 8.5 microns (see Fig 1 
of |Muzerolle et a l. (2009)) at which there is no variability in this 
source and above which the sign of the variability between epochs 
changes sign from that at low wavelengths strongly suggests the in- 
fluence of a global geometric feature such as a thick wall or a warp. 
In short the dominant variability in LRLL 31 object likely comes 
from systematic rather than stochastic effects. 



5 CONCLUSIONS 

For a models of turbulent accretion discs we have derived a pre- 
dicted RPE in the luminosity as function of frequency, focusing 
on the case in which the observational exposure times are small 
compared ti other dynamical time scales. We have discussed the 
implications for interpreting multi-epoch spectral observations of 
spatially unresolved protoplanetary discs. 

The RPE depends on the radial scale of averaging and there 
exists an optimal scale that minimizes the error which also corre- 
sponds to an optimal spectral resolution when the correspondence 
between peak frequency of emission and radius is made. The data 
can be binned to compare with the theory of minimum RPE only 
when the spectral resolution of the telescope exceeds the value 
Ropt which minimizes the RPE. An optimal resolution exists be- 
cause if the spectral resolution is too high then the instrument ends 
up sampling noise and if the spectral resolution is too low, then the 
instrument samples overly coarsely binned regions of the disk. 

The stochastic variability expected between multi-epoch ob- 
servations of discs is an implicit prediction of alpha disc theory and 
can be directly compared to observed variabilities. Complementar- 
ily, the nature of the observed variability can be used to constrain 
whether or not the behavior of the disc is consistent with luminosity 
produced by stochastic or systematic processes. 

The predicted stochastic variability from a discs has several 
distinct characteristics that would signature its prevalence: 

• The stochastic variability is predicted primarily for regions 
where the luminosity is dominated by internal turbulent dissipa- 
tion in the disc. For regions dominated by surface dust-reprocessed 
starlight, the stochastic variability could be expected to be strongly 
reduced. 

• . For regions of the discs dominated by from turbulent dissi- 
pation, the strength of the stochastic variabilty would increase grad- 
ually with decreasing emission frequency (and thus increasing ra- 
dius) down to a critical frequency below which the associated disc 
radii have orbital times exceeding the inter-epoch observation time. 
There the expected variability from the RPE would drop sharply. 

• For inner regions of the discs in transitional objects, where 
emission may be dominated by turbulent dissipation but the discs 
are optically thin, the stochastic variability would be expected to 
be lower than that expected from the same radii for optically thick 
turbulent dissipation dominated discs. 

In applying the basic ideas herein to the specific case of LRLL 
31, we find that systematic variability NOT stochastic variability 
dominates. This is indeed consistent with the conclusion of Muze- 
|rolle et aT] l |2009} . As more multi-epoch observations of broader 
samples of discs at different stages in their lifetimes are obtained, 
the present work may help provide a tool to distinguish stochastic 
processes in discs from systematic dynamical changes and the time 
scales on which these occur. 

Two overall lessons from this analysis are: (1) the precision of 
axisymmetric accretion theories can be quantified and deviations 



r(au) 


a 


^opt 


^ (%) 






0.1 


0.001 


1832 


0.08754 


200 


0.09226 




0.01 


579.3 


0.2781 


200 


0.2816 




0.1 


183.2 


0.8914 


200 


0.8915 


1 


0.001 


1156 


0.1103 


200 


0.1139 




0.01 


365.5 


0.3514 


200 


0.3531 




0.1 


115.6 


1.135 


200 


1.144 


10 


0.001 


729.3 


0.1392 


200 


0.1417 




0.01 


230.6 


0.4449 


200 


0.4451 




0.1 


72.93 


1.454 


200 


1.513 


20 


0.001 


634.9 


0.1493 


200 


0.1515 




0.01 


200.8 


0.4780 


200 


0.4780 




0.1 


63.49 


1.570 


200 


1.661 



Table 2. Optically Thin Case. In column 4, we list the RPE varibilities cal- 
culated using R & f = R op t and and in column 6 the values corresponding 
to the Spitzer telescope low mode R e f = R te i ~ 200. The RPEs are 
computed from Eq. (32| for different r and a and h /r from Eq. l |33} using 
p = 0.6. The optimal resolution is calculated from Eq. {22} as described 
below Eq. l |33| . Note that if R, op t < Rtel data can be binned to compare 
with the a-disc model of minimal RPE. 



r(au) 


a 


^opt 


^-(%) 






0.1 


0.001 


625.8 


0.6935 


200 


0.7111 




0.01 


249.2 


1.116 


200 


1.118 




0.1 


99.19 


1.835 


200 


1.881 


1 


0.001 


557.8 


0.7355 


200 


0.7510 




0.01 


222.1 


1.186 


200 


1.186 




0.1 


88.40 


1.958 


200 


2.030 


10 


0.001 


497.1 


0.7802 


200 


0.7935 




0.01 


197.9 


1.260 


200 


1.260 




0.1 


78.79 


2.091 


200 


2.200 


20 


0.001 


480.2 


0.7942 


200 


0.8069 




0.01 


191.2 


1.284 


200 


1.284 




0.1 


76.10 


2.133 


200 


2.256 



Table 3. same as Table 2 but for optically thick case. The RPEs are com- 
puted from Eqs. (25) and (31) for different r and a and for the optically 
thick version of h/r in Eq. |26| and p = 0.9. 

the alpha-accretion disc theory and observations at a given epoch 
cannot rule out the theory if the deviations between epochs exhibit 
stochastic behavior and fall within the expected RPE of the the- 
ory. The theory may be incomplete but because a mean field theory 
for a turbulent system is intrinsically imprecise, but the impreci- 
sion must be quantified so that the user realizes its limitations. (2) 
Arbitrarily high spectral resolution can lead to misleading compar- 
isons between theory and observation and binning the data may be 
necessary to compare the data with the theory of minimum RPE. 
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